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16 . Abstract • f 

This 'articles consifers the problem of determining, the 
time of occurrence of gamma flares at different points 
in space. This problem arises in the localization of 
sources- of gamma., flares by the triangulation method 
using several spacecraft; (SC). A study is made of different 
methods for determining the time at which the flare reaches 
the . spacecraft and algorithms are .given- for estimating their 
accuracy. The results of processing, the model information 


are also given. 
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Determination of the Time of Occurrence of 
Gamma Flares at Different Points in Space 


B. L. Noval 


1. One of the most interesting results of extra- 
atmospheric astronomy of recent years is the discovery 
of powerful flares of space electromagnetic radiation in 
hard x-ray and gamma-ranges. One very important problem 
of studying the sources of the flares is the determination 
of their position in the heavens (localization). The most 
effective method for solving the problem of localization is 
to record the time at which the flare arrives or, which is 
the same thing, the delay in the arrival of the flare at 
different points in space. This recording is done by sen- 
sors which are carried on the spacecraft (SC) at different 
distances from each other. Thus the determination of the 
celestial coordinates of the source reduces to statistical 
processing of recordings of the flare obtained on different 
SC. 


The localization accuracy depends on the accuracy of 
establishing the position of the SC in space and the accur- 
acy of recording the time of occurrence of the flare by the 
SC. As follows from [1], for real cases of finding the SC 
in space the accuracy of determining their coordinates is 
such that the error in the time of occurrence of the flare 
can be reduced in every case to the time resolution of the 
existing sensors. Taking this fact into account, and also 
the fact that the flare has a finite duration, we must esta- 
blish what we mean by the flare time of occurrence. Let us /4 
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consider the process of recording the flare. Let us assume 
the sequence{_ M j | j 3 * JL,2., ..,M represents the reading of the 
sensor carried on the i-th SC and fixing the calculation (time 
spectrum) with the resolution At . With the very natural 
and unlimited assumptions regarding the nature of the radiation 
and the recording, the sequence J ft . | may be regarded as a 
Poisson process with a certain variable, which is unknown 
previously, with the intensity . If we use Identical 


sensors, then the functions 




for different values of i 


differ only by the displacement of the arguments, and thus we 
may speak of the function A(*) , having in mind all A (•) 

Thus, the problem of determining the delay for the flare to 
reach 2 SC (we shall call them the zero one and the first one) 


is reduced to finding the , such that A (ii*T)='A («V 

X}< 


If the time spectra 


K 


rcoincide, then the prob- 


lem of identifying them, and consequently determining X , 
would be trivial. However, such an agreement does not occur, 
due to the discreet nature of the recording process and the 
random nature of the radiation process. 


The sensors which are designed to study the flares usually 
carry special discharge equipment which establishes when the 
threshold value has been reached in a certain period of time. 

The signal "flare" is thus given. The arrival time of the flare 
may be determined in principle, using this equipment. Let us 
determine the error At in determining the delay when using 
this equipment. We shall confine ourselves to examining only 
the errors due to the discreet recording of the radiation, i.e., 
without considering statistical fluctuations. Using a linear 
approximation in the estimation, we obtain /5 

AX — * 


XTOT3 


or 



A 0>S1> ; , 


A 


(I) 
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) 



where A , K are the parameters determining the linear 
approximation of the function A(-) i T 0 - value of the 
argument of the function M') , corresponding to the beg- 

inning of the interval in which the recording reaches the 
threshold value on the zero SC; - difference in the 

values of the arguments of the function MO , corresponding 
to the beginning of the time intervals in which the record- 
ing reaches the threshold value on the zero and the first 
SC; - quantity, determined by the logic of the dis- 

charge equipment operation. 


We shaii^use K to designate the integer estimate of 
the ratio . Using the formula obtained above (I) and 
the characteristics of the equipment "SNEG-2M" [2], we con- 
clude that K may reach 8. This Is an impermissibly large 
quantity. We therefore turn to the methods based on the 
more complete use of the measurement information. We should 
note that K - the integer estimate of the relationship -A'E' s 
determined by (I), will be used below. 

2. Each of the 1 = lg the 

realization of the random quantity .Jx j , distributed 
according to the Poisson law with a certain parameter. It 
is desirable for us to have the sequence of realizations 
of normal random quantities with unit dispersions. These 
sequences may be obtained by the following normalization /6 


tCi) 


n 


ft) 
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< D 






This normalization is based oxuJt^e estimate of 
the mean square deviations y . For j'V* 

the following interval with ^iven conned 


^ f ° r 
jH / we may give 

& iven conricience coefficient 


\\l - cc^Vn 7 * 
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W , (i) . (i) vlfui' 
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whose magnitude depends on the numb^ Cfc) [33. Therefore 
a study of the difference may be replaced by 

the study of w 


11 1 


n? 


f~w 

0 A 


, * cW'fhf 


Using theL’Hopitalrule, we obtain 

■ U J ffl - £& 


,itn 

Q 


i 2 


This result shp ws the validity of using 


Vnl 


an estimate for a. However, to decrease the scatter 

T J L 

* Jf 1 I 5 we first smooth the estimate 


6 


as 


'6 


for finite values u*j 
taking 5 points for the 'smoothing. Thus we obtain the se- 
quence 


A 


/7 


wh'-;re 


at 






Forming thpsp ^cmennes. we shall look for the minimum 

^lf/J (•) j - 
1 Sv i °^2-K (2) 

The minimum of this expression may be written with respect 

to the index t . The value, which represents the minimum 

0 « 

will be designated by »• * 

0 


Summation in (2) is performed with respect to the index 
S from I to m, where m equals al-K ; K - the integer estimate 
A'F determined by the expression (I). 

At 


Let us examine the accuracy of identifying the spect ra , 
i.e. finding , done in accordance with (I). It is clear 
that when there are no random fluctuations the accuracy would 
be determined by the quantity *t- Actually when finding 
*£* we replace the solution of the continuous problem having 
the form 



by the solution of the discreet problem 
m. 

E(}. x (SAt^t)-4(SAt))- 





for fnnntlons UO’U-l such that there is a ^ for 
whic Under the assumption that At 18 

small with respect to the time for a large change in U-) 


and 


the error of this replacement is At 


/a 


Below we shall estimate the accuracy of identification 
taking into account the fluctuations. Since there is no 
important a priori information regarding the function M-) 
the estimates will be a posteriori. 


3. The first method of accuracy estimation. In accord- 
ance with the given normalization, we will calculate each 
of the-v.^- • ^ , by the. realization of a normally distri- 
buted random quantity c with unit dispersion. Therefore, 

' ' V 


is the realization of the random quantity 


a non-central 


\ 


- square distribution with m ai 


having 
• gre e s of 


freedom and a certain non-central parameter*^ 
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Using the normal approximation for the non-central 
X“ square distribution, * f e formulate the confidence inter- 

V m / ^ , % 

In view of the approximation ilJ -«nirrr iiiA 

M2 (tn+ZM;) 

is a normal quantity with parameters of zero and ^ne^ 

If we use ^ to designate the quantile of qrde,^ - 


val for 


of the normal distribution Cl.e. PiHU-V /or the 
normal quantity £ with parameters 01 zcrv and one), then 

P{~^/ f-i~( n,r W Q |and 1-cC. 

V2(m+ai p J «(l ~ 

Here ^ is tne cumxdence coefficient. 

It is apparent that _ c ^ 

1 * !>»,)• v M i 2iin+4V; * 

Thus we obtain the relatloship for the confidence values of 
the parameter V.* 


and the expression for the upper and lower 

confidence boundaries. A 






+1 — tn 
a 


(3) 


At first glance it would appear that the expression 
under the root sign in ( 3 ) may be negative and consequently 
real confidence boundaries may not exist. 


We shall show that by considering the probable meaning 
of the quantities included in this expression and making 
certain natural assumptions, we can avoid this situation. 


Let us consider the solution of the problem 
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Y+m- CC^Vi(m +2^) '-*• ruin. 


(4) 


under the condition 


V« 


The quantity^ etermines, with the 
selection of the corresponding CO*) , the minimum possible 
value of the random quantity having the non-central distri- 
bution - square with the non-central parameter "Vj . 

The solution of problem (4) is zero at 2C((0^m 3 ’ Usually 
in applications of the theory there are not more than 4 
probabilities Cfc) and therefore the solution is 
for any reasonable number of measurements. 


zero 


Thus, the minimum value of the random quantity having 
a non-central - square distribution with m degrees of 
freedom for an unknown , is determined by the expres- 
sion fT\ ” and ifT" i 

-C(A)V2.m 

It is apparent that if 

h-c«|£m < 

then the expression under the root sign in (4) is not 
negative. 


It follows from (5) that 



/1C 


( 6 ) 


then the inequality is satisfied for all \ ^ everywhere. 


The inequality (7) leads to the condition 
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(n > 8cfo 


and, again considering the probable meaning of cM , we 
obtain 


18) 


If it is not necessary to satisfy the condition C7), then, 
considering that the solution of the problem 


% - awVa?-* T 1 ’ 


equals the minus 


CM 


we reach the condition 




(9) 


rn)Ll% 


or if at m>a , 


Thus, ifHl/145 for an arbitrary 
then there are real confidence boundaries ^ . We shall 

assume that any of these unlimited conditions are satisfied. 


Using the confidence intervals formulated, we obtain 
an estimate of the accuracy for determining . For this 
purpose we use the system of intervals 

and examine , where those 

are included for which the following condition is satisfied 


T. 


* 




( 10 ) 


The estimate of will be determined by the condition that 

't belongs to the set ( ^ . 


The second method of estimating the accuracy. To 
improve the estimate of accuracy of Identification, it is 
necessarjv to consider the stochastic dependence of the quan- 
tities for different values of ^ , for which it is nec- 

essary to Know the concurrent distribution of the quantities . 
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For this purpose we obtain the expression for the character- 
istic function of quadratic form z~X(\X “ here X 15 
normal with M«4f‘Zi4M 


rx-p. 


unit matrix) 


•} 


( 11 ) 


The density function of the quantity f, is determined 
by the expression "X we obtain 
the following for the characteristic function of the random 


quantity 


and 


tf (t) =^)" T i •• \ (xLtfi)exp (-^ (xyfp-fidx, 


( 12 ) 


where i is imaginary unity. 


Ii is apparent that the integrand equals 


Let us transform it to the form 

(x+U) yf (x+li) + v 

for U .V Y we obtain the relationships 

£V/U=i« , _ 

U T WU+V“-^ 


(13) 


/12 


Thus we have 




After transformation to the form (13) we may readily see 


9 


that the Integrand In (12) represents the product of the 
distribution density of the normal quantity with the math- 
ematical expectation minus and the covariational matrix 

V* 1 ‘ ,nd (yjj4 • Thue 

Kfl - |E-2 itftf V 4 ^ 

or 

(1 „ 

Using expression (14) we may readily find the first 
two moments of the random quantity j? 

it ~ x f '(°) = Sp ft > 

& 4 -k^ o M/ , V + V < 0 ,+ U5> 

a S 

$$ are the diagonal elements of the matrix A . 


T hUS 

Let us set 



where ■'j 
its element 


/ 1 3 

. is the dimensionality matrix * N| , such that 
^ ^ equals unity at £ ^ ^ rn t and the re- 


J 
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maining elements equal zero 


/. _-!/o r 

Then satisfies the conditions ( 11 ) and 2 i . 

ifference £: and £* t vie ^he quadratic form 5* 

.1 A zf **■ 1 \ • 


in view of the assump- 


The difference *», . 

, «J y o jl 

On the other hand, 2 L ’ and 
tions made have § non-central \ - square distribution. 

Therefore <&,-*/> equalE ■ , -'X''V) ’ and i (i, 
according to ( 15 ), equals 1 ^ ^1 


*\\ 


K 


c’Wj-BfV -i 


Let us assume 


... . 

\if (ft r rt r Vr ? fi M i ' ) ?; 


(16) 


W 


We shall assume that a normal approximation is valid for 
the difference of the quadratic forms Included in the denomina- 
tor, Then expression ( 16 ) represents the realization of the 
normal random quantity with zero mathematical expectation and 
unit dispersion. 

Since is not known, we make an estimate for^ . 

We may readily obtain the following from (15) 


T ^ j ' / 

=i.‘ ( f)j as the estimate for 

•'%-A* r- We shall assume that replacement in (15) 

.j _ _ ..'IV _i j i l 1. 1 ■» .. m - t j 


And we use-1 


does not disturb the normal approximation. 

. jLZHi* - v ‘ V c 7# 

Let us set ^ 

Using the statistics *iU). , we follow the procedure for es- 
timating the accuracy, based on verifying the statistical 


li 



v * —V s a 

hypothesis <J . To solve the problem of verifying 

this hypothesis, it is sufficient to use the results of the 
most valid non-biased criteria for unilateral alternatives /14 
for the exponential sets £ 33 * Our problem is a very partic- 
ular case of this more general problem. Using these results, 
we obtain the following procedure. 


I 


1- 


Let us assume is the quantile of order of the 

normal distribution; ^ - error of the second kind; /j 

, as in the previous section. 


Then, if 


t j(a-) 


f.HA ('nr>T>p 


i, 


for Wnicn the 


J: 




must be verified. 


(17) 


pnr>T^>sponding interval is included in the union 

It ^is apparent that only those values of ^ 


corresponding to them are included in 


The error is determined the same way as was done 

3n the previous section, by the condition of belonging to 
the set jt) . We should note that both the set p> , 

P** may be comprised of non-intersecting seg- 


and the set 


merits. This case may be interpreted as the existence of sev- 
eral solutions of the problem for determining 


5. Computational algorithms and computer programs were 
formulated using the procedures given above. 


In formulating the algorithms, we considered the possible 
singularities having the form of measurement information, 
particularly the occurrence of inc' rrect measurements. Mea- 
sures making it possible to process this information were 
taken into account in the algorithm. 
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The programs were used to process the results of Earth- 
based tests of the "SNEG-2M" equipment. This equipment they 
used for the Soviet-French experiment on studying the flares /15 
of gamma radiation in space [2], The total number of measure- 
ments M by the equipment equalled 1024. This was sufficient 
for the use of the approximations employed in the procedures 
mentioned.. The measurement information was obtained by using 
a flare model. 

We shall give the results of processing the measurement 
information. The minimum value of ^ was obtained when the 
indexes in the first and zero time spectra coincided. This 
pointed to the fact that in this case the signal "flare” 
was processed very accurately. 

When estimating the accuracy by the first method, the 
confidence level was assumed to equal 0.9. This moderate 
value was selected due to the large rigidity of the solution 
rule. The use of the procedure given in this method resulted 
in 9 intervals of J , for which the condition (10) was 
satisfied. The was the interval of 

the quantity 13 At , and the 1 * possible deviations from the 
most probable solution were determined by the quantities 
4c * ^4 and 

f w 

This accuracy of determining £ was Inadequate and there- 
fore the second method of accuracy estimation was used. Based 
on the same considerations as in the use of the first method, 
the error of the second type was chosen to equal 0.1. The 
application of this procedure resulted, in the formation of 
two mixed intervals 3^ and the set the interval 

of the quantity 3st . The possible deviations from the 
most probable solution/^ vaai, 4 1 . 
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